^AoB Plants 



The open-access 
journal for plant sciences 



r 



Research Article 



Spatial and temporal functional changes in alpine summit 
vegetation are driven by increases in shrubs and graminoids 

Susanna Venn 1 - 2 *, Catherine Pickering 1 and Ken Green 3 

1 School of Environment, Environmental Futures Centre, Griffith University, Gold Coast, QLD 4222, Australia 

2 Department of Botany, Research Centre for Applied Alpine Ecology, La Trobe University, Bundoora, VIC 3086, Australia 

3 New South Wales National Parks and Wildlife Service, Snowy Mountains Region, PO Box 2228, Jindabyne, NSW 2627, Australia 

Received: 14 October 2013; Accepted: 24 January 2014; Published: 21 February 2014 

Citation: Venn S, Pickering C, Green K. 2014. Spatial and temporal functional changes in alpine summit vegetation are driven by 
increases in shrubs and graminoids. AoB PLANTS 6: plu008; doi:10.1093/aobpla/plu008 

Abstract. Classical approaches to investigating temporal and spatial changes in community composition offer only 
partial insight into the ecology that drives species distribution, community patterns and processes, whereas a func- 
tional approach can help to determine many of the underlying mechanisms that drive such patterns. Here, we aim to 
bring these two approaches together to understand such drivers, using an elevation gradient of sites, a repeat species 
survey and species functional traits. We used data from a repeat vegetation survey on five alpine summits and mea- 
sured plant height, leaf area, leaf dry matter content and specific leaf area (SLA) for every species recorded in the 
surveys. We combined species abundances with trait values to produce a community trait-weighted mean (CTWM) 
for each trait, and then combined survey results with the CTWMs. Across the gradient of summits, more favourable 
conditions for plant growth (warmer, longer growing season) occurred at the lower elevations. Vegetation composition 
changes between 2004 and 2011 (according to non-metric multi-dimensional scaling ordination) were strongly af- 
fected by the high and increasing abundance of species with high SLA at high elevations. Species life-form categories 
strongly affected compositional changes and functional composition, with increasing dominance of tall shrubs and 
graminoids at the lower-elevation summits, and an overall increase in graminoids across the gradient. The CTWM 
for plant height and leaf dry matter content significantly decreased with elevation, whereas for leaf area and SLA it 
significantly increased. The significant relationships between CTWM and elevation may suggest specific ecological 
processes, namely plant competition and local productivity, influencing vegetation preferentially across the elevation 
gradient, with the dominance of shrubs and graminoids driving the patterns in the CTWMs. 

Keywords: Alpine vegetation; community composition; functional composition; functional traits; GLORIA; 
Snowy Mountains. 



Introduction 

In the context of environmental change, linking quantita- 
tive measures of plant species physical characteristics, 
for example their functional traits or life form, with spe- 
cies distributions and local environmental factors can 
reveal the processes that drive patterns in vegetation 
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composition (Diaz and Cabido 1997; Tilman et al. 2001; 
Lavorel and Gamier 2002; McGill et al. 2006). The inter- 
action among functional groups of species with climatic 
changes undoubtedly affects community composition 
and underlying ecosystem functioning (Hooper and 
Vitousek 1997; Westoby and Wright 2006; Spasojevic 
and Suding 2012). Hence, measures of functional traits, 
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functional diversity or functional composition, in combin- 
ation with measures of species composition and diversity, 
may provide greater insights into the assembly processes 
that drive community composition and for gauging 
ecosystem stability in changing environments (Tilman 
et al. 1997). 

Using environmental gradients is a simple way to 
examine natural variation in vegetation and community 
responses to environmental changes (McGill et al. 2006). 
Intraspecific variation and species turnover can cause 
shifts in trait values differentially across gradients, and 
thus functional diversity across environmental gradients 
is a function of biotic and abiotic interactions (Venn 
et al. 2011; Schbb et al. 2012; Spasojevic and Suding 
2012), disturbance regimes (Choler 2005), population dy- 
namics (Pollock et al. 2012) and biogeochemical cycling 
(Mason et al. 2012) at local, regional and biome scales 
(Reich et al. 1997). Mountainous regions are therefore 
ideal for testing the role of plant functional traits across 
environmental gradients in determining community 
composition and shifts therein; over relatively short spa- 
tial scales, rapid changes in elevation interact with local 
topography to create steep gradients in temperature 
and precipitation (Korner 1999). Thus, the variation in 
species abundance across an environmental or elevation 
gradient will impact on community composition and 
ecosystem functioning (Tilman et al. 1997), as different 
processes are affected by different species via changes 
in the representation of species functional traits (Chapin 
et al. 1996). The interrelation between a species life form 
and its functional traits is exemplified across elevation 
gradients in alpine regions; for example, shrubs (taller, 
long-lived, woody plants) dominate in favourable envir- 
onments. Through their morphological and physiological 
traits, shrub species can modify a wide range of ecosys- 
tem processes, including alteration of local snow depths 
and associated hydrological dynamics, nutrient exchange 
and associated net carbon balance (Myers-Smith et al. 
2011). In addition, shrub species are often taller than 
neighbouring forbs and they can be competitively super- 
ior, forming dense thickets with closed canopies. In- 
creases in shrub cover and height can also potentially 
restrict the growth of other plant species by limiting 
light availability (Myers-Smith et al. 2011). Graminoids, 
namely grasses and sedges, however, may compete 
with shrubs in the more favourable environments, recruit 
within the canopy of senescing shrubs (Williams and 
Ashton 1988) and are well adapted to a range of environ- 
mental conditions; many C3 grasses and sedges can also 
rapidly increase in abundance after sufficient rainfall and 
improved abiotic conditions (Jarrad et al. 2009). 

Here, we use community trait-weighted means 
(CTWMs) to explore the driving mechanisms of temporal 



change in vegetation composition across an elevation 
gradient of alpine summits. We used four easily mea- 
sured morphological traits: plant height at maturity, leaf 
area, leaf dry matter content (LDMC) and specific leaf 
area (SLA); some of these are interrelated, but all repre- 
sent important dimensions of functional and strategic 
variation among plant species (Westoby 1998; Weiher 
et al. 1999) and are responsible for some of the most strik- 
ing and important patterns in species distributions in the 
field (Westoby et al. 2002; Choler 2005). In addition, 
large-scale comparisons among biomes provide evidence 
that these traits may be viewed as relevant functional 
markers suitable for predicting species performance 
along gradients (Reich et al. 1997). While we did not 
measure community assembly processes directly, these 
chosen traits are important for various assembly pro- 
cesses and population dynamics such as competition, 
facilitation, productivity, and stress tolerance and longev- 
ity (Grime 1977; Westoby 1998; Cornelissen et al. 2003; 
Schbb et al. 2012; Michalet et al. 2014). We can therefore 
infer the mechanisms behind changes in community 
composition based on life forms, individual traits and 
patterns in the CTWMs. The elevation gradient of sites 
represents a strong gradient of temperatures, snowmelt 
date and growing season length, and therefore repre- 
sents an important gradient across which community 
composition and functional composition are expected 
to vary preferentially, according to position along the 
gradient. Specifically we ask: what is the spatial and 
temporal variation in species composition across the 
gradient of sites; and where are various combinations 
of traits and life forms most prevalent? We then discuss 
how these composition patterns and functional traits 
interact with species' ecology preferentially across the 
gradient of sites. 

Methods 

Study sites 

We used five alpine summits in this study that represent 
the Australian contribution to the Global Observation Re- 
search Initiative in Alpine Environments (GLORIA): an on- 
going, empirical study with the specific aim of detecting 
alpine vegetation change on summits in relation to cli- 
mate change. In January 2011, we re-surveyed the sum- 
mits that were originally marked and surveyed in January 
2004 that exist along a continuous ridge from close to the 
valley floor to the summit of Mt Clarke (Pickering et al. 
2008; Pickering and Green 2009) (Fig. 1). The summits 
cover an elevation range of 301 m from the lowest at 
1813 m (Clarke 5) through to the highest at 2114 m 
(Clarke 1) (Table 1) and cover a horizontal distance of 
1600 m. The summits were selected for long-term 
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Figure 1. Location of the Snowy Mountains in Australia and study site locations at Mt Clarke 1-5 (CL1, CL2, CL3, CL4, CL5) representing the five 
summits. 



monitoring under the GLORIA sampling protocols (Pauli 
et al. 2004), as they experience similar effects of exposure 
and differences in climate that are most likely due to the 
elevation gradient. They are all relatively flat, rather than 
cone-shaped peaks, and the vegetation is characteristic 
of nearby summits. The soils are around 350 + 110 mm 
in depth (K. Green, unpubl. res.), well-formed alpine 
humus soils (Costin 1954). The highest summits are domi- 
nated by tall alpine herbfield (a mixture of forbs and gra- 
minoids), whereas the lower summits are characterized 
by shrubs. As a result of the continuous, mostly perennial 
vegetation cover, biomass is high compared with some 
other alpine regions (Costin 1954). There are some rock 
outcrops, but these are not a defining feature of the 
summits. Disturbance is minimal as cattle grazing ceased 
>60 years ago, the historical stock travelling route 
avoided these summits and there are few native and 
no exotic burrowing mammals at these elevations. 
No walking tracks cross the summits, resulting in low 
visitation rates. 

Variation in abiotic variables across the gradient 
of summits 

The variation in soil temperatures, nutrients and minerals 
across the elevation gradient has been described previ- 
ously for the summits (Pickering et al. 2008; Pickering 
and Green 2009; Venn et al. 2012a). Temperature record- 
ing began at the summits in January 2002 using 



temperature loggers (Tinytag Plus— Gemini Data Loggers, 
Chichester, UK) buried 10 cm below the ground surface. 
From January 2004 to the present, loggers have been in 
position at the 'corner' of each aspect of each summit. 
Temperatures were recorded every 2 h. Temperature 
data from 2004 to 2011 were used to calculate annual va- 
lues of absolute minimum soil temperature, annual daily 
mean soil temperature, absolute maximum soil tempera- 
ture, temperature sums (>5 °C), growing degree days 
and the length of the growing season across the years 
sampled. Several climate parameters were derived from 
these data and used in previously published analyses 
(Pickering and Green 2009; Venn et al. 2012a). 

Precipitation data were collected from 2003 to 2011 
from a Bureau of Meteorology automated weather sta- 
tion about 8 km to the south at Thredbo (1957 m), and 
from Pengilley Bog (1730 m, during the growing season 
only), 13 km to the north-east. 

Soil samples were collected in 2003 at each summit, 
on each aspect close to the 1-m 2 quadrat clusters, giv- 
ing a total of four individual soil samples per summit. 
A minimum of 500 g was collected in a single sample, 
taking the top 10 cm of soil through a 75-mm core. 
Samples were air-dried, sieved through a 2-mm mesh 
and analysed by the CSIRO Division of Soil and Water, 
Canberra, Australia. Samples were analysed for 17 
minerals/nutrients (see Pickering and Green 2009 for 
full details). 
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Table 1. Summit details including locations, elevation (metres above sea level (m.a.s.l.)) and area (m 2 ); the mean species richness and overlapping vegetation cover for each of the upper and 
lower SASs divided into three life-form groups (graminoids, forbs and shrubs) for 2011 and the change in mean species richness since 2004 (increase + or decrease -) at each summit; and 
the CTWM for four traits (height, leaf area, LDMC and SLA) for each of the sampling units for 2011 and the change in overlapping cover (%) since 2004 at each summit. 
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Vegetation sampling 

The top section of each summit was divided into eight 
summit area sections (SASs), four covering the area 
down to 5 m below the summit, the 5-m isoline, for 
each of the four cardinal compass bearings (hereafter re- 
ferred to as the upper or 5-m SAS), and another four cov- 
ering the four compass bearings down to the 10-m isoline 
(hereafter referred to as the lower or 10-m SAS) (Pauli 
et al. 2011). Where the summit was exceptionally flat, 
the upper sampling area extended 50 m from the summit 
and the lower area extended 100 m. In each of the eight 
SASs the percentage cover of each vascular plant species 
was estimated using a random step-pointing method 
(200 points per area). In this method, vegetation is 
sampled at randomly generated points across the site, 
and hits are summed for each species in order to provide 
a measure of species abundance. Any species not 
sampled by the step-point method were visually assessed 
for cover (Pauli et al. 2011). 

Sampling in 2011 was conducted 'blind' without refer- 
ring to the 2004 data, and was performed by the same 
people to ensure consistency in the data over the days 
of fieldwork (Vittoz et al. 2010). Postsampling, a rigorous 
species identification checking procedure, was used to 



ensure that changes in species between 2004 and 2011 
were accurate. All species names follow Costin et al. 
(2000) to be consistent with the initial 2004 survey. 

We selected four important, character-based plant 
traits (plant height, leaf area, LDMC and SLA) (Table 2) 
to use in the functional composition analyses. Trait 
data, including destructive sampling, were not directly 
measured from plants within the study sites; rather all 
traits were collected from adjacent areas <200 m from 
the study sites, from within the Kosciuszko National 
Park alpine area in growing seasons of 2010, 2011 and 
2012, using 10 adult individuals for each measurement. 
We chose these traits based on the individual plant and 
ecosystem functions in which they are involved (Leps 
et al. 2006; Petchey and Gaston 2006), in the context of 
alpine plant community dynamics (Table 2). 

Data analysis 

Vegetation composition. In order to determine shifts in 
species composition and abundance between 2004 and 
2011 across the gradient of sites, we used non-metric 
multi-dimensional scaling (NMDS) and ordination. 
Species abundance data from the upper and lower 
SASs were log + 1 transformed to down-weight very 



Table 2. Plant traits measured for all species recorded in the upper and lower SASs. Traits were measured in the field and laboratory based on 
protocols outlined by Cornelissen et al. (2003). 

Trait (unit) Description Functional indicator 

Plant height (mm) Shortest distance between the upper boundary of the A measure of species overall competitive ability at plant 
main photosynthetic material (usually the canopy) and maturity. Species that are relatively taller will be more 
ground level competitive, usually for light. Indirect measurement for 

biomass, lateral spread, rooting depth and leaf size 

Leaf area (mm 2 ) One-sided projected surface area of an average leaf A measure of stress tolerance. Small leaves tend to be 

favoured under heat stress, cold stress, drought stress and 
high-radiation stress. Within a climate zone, leaf size tends to 
increase with plant height and soil nutrients, but decreases 
with disturbance. Larger leaves are expected in more 
productive landscapes 

LDMC (mg g -1 ) The ratio of the oven-dry mass of a leaf to the fresh Indirectly represents the mean density of leaf tissue, relates 

weight of the leaf to the inverse of SLA. Low LDMC can indicate fast resource 

acquisition. Leaves tend to be more resistant to physical 
stress such as wind and hail. Species with low LDMC tend to 
be associated with highly disturbed environments and high 
productivity 

SLA (mm 2 mg™ 1 ) Ratio of one-sided area of a fresh leaf to its oven-dry Low values correspond to relatively high investments in 

mass defences to harsh conditions including long life spans and 

structural adaptations. Reflects the expected return on 
previously captured resources such as light and nutrients. 
A good positive correlate of potential growth rate 
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common species, and dissimilarities between all pairs of 
aspects from the upper and lower SASs were calculated 
using the Bray- Curtis dissimilarity coefficient (after 999 
iterations), previously noted as a robust measure in 
recovering ecological distance over a range of models 
and stochastic variations in the data (Bray and Curtis 
1957; Quinn and Keough 2003). How good the 
dissimilarity matrix is was determined by the stress 
value (Kruskal 1964; McCune and Mefford 1999). Stress 
values <0.2 are recommended (Clarke 1993) as values 
above this threshold may mislead interpretations. 

The difference in vegetation composition between 
2004 and 2011 was analysed by an analysis of similarities 
(ANOSIM) procedure (Clarke 1993; Quinn and Keough 
2003). This procedure is analogous to an analysis of vari- 
ance comparing between-group and within-group vari- 
ation. The ANOSIM procedure tests the null hypothesis 
that there are no differences between a priori defined 
groups (in this case years), or that the average rank of 
dissimilarities between all possible pairs of objects in dif- 
ferent groups is the same as the average rank of dissimi- 
larities between pairs of objects in the same groups 
(Quinn and Keough 2003). All ANOSIM procedures for 
abundance data from the three spatial scales used per- 
mutation/randomization methods on a similarity matrix 
to allocate objects randomly to groups and then generate 
the distribution of R under the null hypothesis that all ran- 
dom allocations are equally likely (Clarke 1993; Quinn and 
Keough 2003). The R distribution is scaled between pairs 
of objects in the same group with values between - 1 and 1. 
Differences between groups would be suggested by 
R values greater than zero where objects are more dis- 
similar between groups than within groups. R values of 
zero indicate that the null hypothesis is true. R = 1 indi- 
cates that all samples within groups are more similar to 
each other than to those in different groups (Quinn and 
Keough 2003). 

The dominant character species, those that are useful 
in discriminating between years at each of the spatial 
scales, were explored using the similarity percentages 
(SIMPER) procedure (Clarke 1993), performed using the 
PRIMER 6 package (Plymouth Routines in Multivariate 
Ecological Research 5.1.2. 2010). This procedure utilizes 
the similarity and dissimilarity indices between all pairs 
of samples to identify typical species within a site/year, 
as well as important species that distinguish between 
years. 

Community trait-weighted means. In order to illustrate 
functional changes between 2004 and 2011 across the 
gradient of summits, we calculated the CTWM for the 
assemblages of plants recorded in the upper and lower 
SASs at every summit, for each of the four traits. The 



CTWMs were calculated by the method proposed by 
Mason et al. (2003), interpreted by Mason et al. (2005) 
and modified by Leps et al. (2006). This method allowed 
us to calculate the relative abundance of individual trait 
responses to the gradient of sites and over time, 
weighted by the absolute individual species abundances, 
rather than provide an aggregate index based on 
multi-trait space (Villeger et al. 2008). We utilized the 
software of Leps and de Bello (2008) to assist in scaling 
our trait and abundance values for the data from the 
three spatial scales and to calculate the CTWM values. 

Temporal change across the summits in CTWM was 
measured as the per cent difference between mean 
values between 2004 and 2011, and by examining the 
overlap of 95 % confidence intervals from a normal 
distribution of the data across the three spatial scales 
(Cumming and Finch 2005). Combining all sites, the dif- 
ference between the 2004 and 2011 values of functional 
composition for each trait was compared with paired 
t-tests using the quadrants of each site (aspects) as repli- 
cates. The relationships between site elevation and the 
abundance-weighted trait means were assessed with 
simple linear regression. Additionally, the relationships 
between the vegetation composition across sites and 
trends in the CTWMs at those sites were illustrated by 
overlaying vectors on each of the two-dimensional ordi- 
nations, for both the upper and lower SASs, in order to il- 
lustrate which traits are more abundant and display any 
trends across the ordination, as well as their direction and 
influence in relation to the compositional floristic data. 
Only the community-weighted trait means were used 
as vectors, and only those with correlations to each ordin- 
ation where R 2 > 0.5 are displayed. 

Ordinations, ANOSIM and SIMPER routines were per- 
formed using the PRIMER 6 package (Plymouth Routines 
in Multivariate Ecological Research 2010). SYSTAT ver. 10 
(Copyright SPSS Inc., 2000) was used for all other statistical 
analyses. 

Results 

Variation in abiotic variables across the elevation 
gradient 

Across the gradient of sites, soil temperature and growing 
season length have been shown in previous analyses to 
vary predictably with elevation. Repeatedly, more favour- 
able conditions for plant growth have been found at the 
lower-elevation summits, with cooler, shorter growing 
seasons with (often) less available nutrients at the higher 
summits (Pickering and Green 2009; Venn et al. 2012a). 
As previously described in Venn et al. (2012a), there 
were no consistent summit-specific trends (increases 
or decreases) in climatic variables across the summits 
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between 2004 and 2011, although relationships between 
summit elevation and climate between 2004 and 2011 
remained consistent; minimum temperatures (linear 
regression: R 2 = 0.28, P = 0.016), mean temperatures 
(linear regression: R 2 = 0.44, P = 0.005) and number of 
growing season (snow-free) days (linear regression: 
R 2 = 0.23, P = 0.032) all significantly decreased with in- 
creasing elevation but not maximum temperatures 
(using the data from four temperature loggers at each 
site as replicates) (Venn et al. 2012a). 

Although precipitation was not specifically measured 
at each summit, data from the nearby rain gauges re- 
vealed substantial increases in annual and growing sea- 
son precipitation over the 2010/2011 growing season 
(October to April) in comparison with previous years in 
which the region experienced low rainfall conditions 
for almost a decade (Venn et al. 2012a, b). At Thredbo, 
there was a 30 % increase in annual precipitation in 
2011 compared with 2004 (mean precipitation between 
2004 and 2009 was 1136 mm, whereas in 2011 it was 
1647 mm). Growing season precipitation at Pengilley 
Bog was about 50% higher in 2011 compared with 
previous years (mean growing season precipitation be- 
tween 2003/2004 and 2009/2010 was 575 mm, whereas 
in 2010/2011 it was 1182 mm) (Venn et al. 2012a). 

Soil nutrient analyses conducted in 2003 revealed few 
trends across the elevational gradient of sites (Pickering 
and Green 2009). Importantly, however, organic carbon 
(%), available nitrogen (as ammonium) and aluminium 
(mg kg -1 ) all decreased with increasing elevation. 

Temporal change in vegetation composition 

Overall, species richness between 2004 and 2011 tended 
to increase at all sites (Table 1), with the largest mean in- 
crease occurring in the lower mountain SAS at Clarke 4 
(+11.5 species). Vegetation composition changes be- 
tween 2004 and 2011 across the gradient of sites were 
most pronounced among graminoids and forbs at the 
lowest site, Clarke 5, with an increase in forb abundance 
of up to 350% in the lower mountain SAS (Table 1). Not- 
able increases in graminoids were partly made up of the 
presence of new species in 2011 in many SASs across 
all summits, namely Agrostis sp., Deyeuxia crassilisica, 
Australopyrum velutinum (Poaceae), Carex breviculmis 
(Cyperaceae) and Luzula alpestris (Juncaceae). 

Changes in abundance across life forms were generally 
within the same order of magnitude between the upper 
and lower SASs on each summit. Shrub abundance in- 
creased substantially in the lower SASs at the two highest 
sites, Clarke 1 and 2 (Table 1). Simple linear regression re- 
vealed few discernible patterns in the change in abun- 
dance of particular life forms between years across the 
gradient of sites, except for within graminoids in the 



upper and lower SASs, with per cent change decreasing 
significantly with site elevation (upper SAS: R 2 = 0.92, 
P = 0.008; lower SAS: R 2 = 0.99, P = <0.001). 

Compositional changes between 2004 and 2011, illu- 
strated by the ordination diagrams, point to similar direc- 
tional changes across all summits and sampling areas. 
In general, the vegetation composition in 2011 has chan- 
ged in a consistent manner within and among sites in 
ordination space (Fig. 2), tending towards groupings of 
species with similarly high SLA. This also indicates that 
the turnover and abundances of key species have been 
similar across sites. This is most pronounced in the dia- 
gram illustrating the combined summit data for the 
lower SASs, indicating that the vegetation composition 
between 2004 and 2011 changed similarly across sum- 
mits (Fig. 2B). The ANOSIM results point to significant dif- 
ferences in the vegetation composition between sample 
times, using aspects as replicates across sites (Table 3) 
within quadrats, the upper and lower summit areas. 

The vectors that represent the functional composition 
in the community-weighted trait means, overlaid on 
the ordination diagrams, point clearly to the higher- 
elevation summits having high abundances of species 
with higher SLA ratios and larger leaf areas (Fig. 2). Con- 
versely, the lower-elevation summits appear to have 
proportionally higher abundances of species that are tal- 
ler and have larger LDMC values. Graminoid and shrub 
abundance tended to contribute to the majority of the 
differences between 2004 and 2011, revealed through 
the SIMPER analyses, using aspects on summits as repli- 
cates and with the data combined for each summit 
(Table 4). The abundances of the top 10 most influential 
and discerning species generally increased between 
sample times (Table 4), leading to the dissimilarity in 
vegetation composition between 2004 and 2011 as 
identified in the ANOSIM analyses. At the lower and 
higher SASs across sites, shrub species were identified 
as the most influential in determining differences be- 
tween years (at least 5/10 of the top species were 
shrubs, using aspects as replicates and with data com- 
bined for each summit respectively) (Table 4). Graminoid 
species were the second most influential group leading 
to dissimilarity between years. 

Spatial and temporal patterns in functional traits 
and CTWMs 

The most substantial differences in the CTWM occurred at 
the lowest summit, Clarke 5, of up to 110% for leaf area 
and SLA traits between years. At the higher-elevation 
summits, the CTWM for leaf area and SLA traits increased 
by up to 26% in 2011, with the largest increases within 
the lower SASs (Table 1). The CTWM for SLA increased sig- 
nificantly between 2004 and 2011 in both the upper and 
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Figure 2. Non-metric multi-dimensional scaling ordinations of the (A) upper and (B) lower SASs at two levels of replication, within the five Mt 
Clarke summits using aspects, north, east, south, west, as replicates, and whole summits, to illustrate differences in vegetation composition 
(abundance) between 2004 (closed shapes) and 2011 (open shapes). Sites (shapes) close to each other in ordination space indicate high levels 
of similarity in species composition. Vector lines represent the CTWM values for the plant traits (height, leaf area, LDMC, SLA) with a Pearson's 
correlation of >50 % with the vegetation data for each year and sampling unit. Vector lines indicate the degree to which those functional traits 
affect the groupings (similarity) of sites in ordination space. 



Table 3. Results from the ANOSIM procedure across summits at the 
three sampling units and two levels of replication within sites 
(summits), comparing the similarity in vegetation data between 
2004 and 2011. 



Sampling unit/replication 


n 


Global R 


P 


Upper SAS/summit aspects 


10 


0.04 


0.01 


Upper SAS/summits 


5 


0.09 


0.19 


Lower SAS/summit aspects 


10 


0.07 


0.05 


Lower SAS/summits 


5 


0.06 


0.25 



lower SASs of Clarke 5 (1813 m) and Clarke 3 (1992 m) 
(Fig. 3). The CTWM for leaf area was significantly greater 
at Clarke 5 in the upper SASs between 2004 and 2011 and 
the CTWM for plant height decreased significantly 
between 2004 and 2011 at Clarke 5 in the upper and 
lower SASs (Fig. 3). 

Paired t-tests between the CTWM for each trait, com- 
bining all sites and using summit aspects as replicates 



(n = 20 for each test), revealed several significant differ- 
ences between years; on the upper summit areas the 
CTWM for plant height and LDMC were significantly 
lower in 2011 (P= 0.004 and P= 0.003, respectively), 
whereas the CTWM for SLA was significantly higher 
(P < 0.001). On the lower summit areas, the significant 
differences trended in the same direction (plant height, 
P = 0.019; LDMC, P = 0.011; SLA, P < 0.001). 

Across the gradient of summits there were several sig- 
nificant linear relationships with good predictive power 
between the CTWM value for each of the measured traits 
and site elevation. Notably, plant height (decrease with 
elevation), leaf area (increase with elevation), LDMC (de- 
crease with elevation) and SLA (increase with elevation) 
were all significantly related to elevation in 2004, where- 
as only plant height and leaf area were significantly re- 
lated to elevation in 2011 (Fig. 3). Considering species 
preferentially organize themselves across the gradient 
in relation to environmental variables (Pickering and 
Green 2009; Venn et al. 2012a), we are confident that 
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any intraspecific trait variation will be overridden by the 
strong patterns in species distributions and life forms at 
these sites. 

The relationships between site elevation and the CTWM 
values for each trait (Fig. 3) demonstrate how the CTWM 
of each site is determined floristically (by abundance) 
across the gradient of sites. Additionally, the significant 
increases in the CTWM for SLA between 2004 and 2011 
can be attributed to the changes in mean abundance of 
key shrub species (Table 4) with small, tough leaves such 
as Kunzea muelleri which increased in abundance. The 
significant decreases in the CTWM for plant height at 
Clarke 5 (Fig. 3), however, can mostly be attributed to 
the relative increases of the short/prostrate shrub Eparcis 
microphylla (Ericaceae) on this summit. 

Discussion 

Linking changes in vegetation composition 
with functional composition 

Our exploratory analyses have demonstrated how a 
short-term change in vegetation composition affects 
the CTWM and how it varies preferentially across the 
gradient of summits. Not surprisingly, the patterns in 
life forms across sites strongly reinforce the patterns in 
the CTWM, as do broad functional groups of species 
across other environmental gradients (Reich et al. 1997; 
Choler 2005). 

In general, with increasing elevation within both the 
lower and upper SASs, shrub species appear less abun- 
dant and forb species more so, whereas the proportion 
of graminoid species tends to increase in abundance 
with elevation, a trend similarly reported in other local 
studies (Venn and Morgan 2005; Venn 2007). Across the 
gradient, the CTWM for plant height decreases signifi- 
cantly with elevation as the proportion of tall shrub spe- 
cies (particularly Phebalium ovatifolium) decreases with 
longer snow seasons and lower temperatures at higher 
elevations. Similarly, the CTWM for the leaf traits corre- 
sponds with the higher proportion of forbs and grami- 
noids at the higher elevations; SLA and leaf area 
increase with elevation and LDMC decreases, reflecting 
the higher abundances of smaller, compact forb and gra- 
minoid species (particularly the sedges Carex hebes, 
C. breviculmis) at the higher elevations, as is found in simi- 
lar studies (Kudo et al. 1999; Choler 2005; Venn et al. 
2012a). These trends were also apparent at the scale of 
SAS within mountains, where directional differences in 
life form and the CTWM exist between the upper and 
lower summit areas. 

Species with small, thick leaves and low SLA are key 
features of high, wind-exposed sites (Tranquillini 1964; 
Kudo et al. 1999; Choler 2005), as is the high abundance 



of short-statured species as found at these sites, but also 
many of the shrub species dominant at the lower eleva- 
tions, which is typical of many Australian shrubs (Groves 
1994). This is consistent with studies relating canopy 
architecture and plant height to wind in exposed sites 
(Caldwell et al. 1974; Smith et al. 1995; Korner 1999) or 
possibly phylogenetic conservatism of habitat (Cornwell 
and Ackerly 2009). In the context of this study, high SLA 
and leaf areas at the higher, supposedly more stressful 
sites appear counterintuitive and are inconsistent with 
similar studies of leaf traits across growing season and 
snowmelt gradients (Kudo et al. 1999; Spasojevic and 
Suding 2012). In the Snowy Mountains, fewer shrubs are 
present at higher elevations and forbs at high elevations 
are generally selected to have smaller leaves, although 
their leaves are generally fleshy and larger than those 
of the shrubs, hence the patterns in SLA and leaf area 
across the gradient. However, at the highest wind-swept 
areas of feldmark, several tough, dwarf, prostrate shrub 
species are present (such as Epacris and Chionohebe 
spp.) in an environment where long-lived species are se- 
lected for. Large-scale interspecific comparisons have 
strengthened the idea that SLA is related negatively 
with leaf life span and positively with relative growth 
rate (Reich et al. 1999; Wright et al. 2004), a potentially 
advantageous strategy in harsh climates with short grow- 
ing seasons. However, a longer snow season at the higher 
summits in the Snowy Mountains may actually provide 
more favourable conditions in early spring when plants 
are protected from early-season frosts and winds, where- 
as they would be exposed earlier at the lower summits 
(Inouye et al. 2002; Venn et al. 2012b), which may have 
caused these unexpected patterns in the CTWM for SLA 
across the gradient. This is certainly the case around 
the edges of late-lying snowpatches in the Snowy 
Mountains, which melt out earlier than the centre of the 
snowpatch (Venn et al. 2011), and where tough, high 
fibrous (low SLA) species are more abundant in the 
early-exposed alpine areas (Choler 2005). Conversely, 
low LDMC and large leaf areas at higher elevations may 
be indicative of fast resource acquisition and relatively 
high productivity, even in areas where growing seasons 
are short and unpredictable (Kudo et al. 1999). 

Only a few studies have documented temporal func- 
tional diversity or functional composition change in 
plant communities, with mixed results; but see the review 
in Diaz and Cabido (2001). Here, significant differences in 
the CTWM of several traits between years, but some non- 
significant differences in terms of species composition 
(from ANOSIM analyses) point to species abundance 
change, rather than high turnover or increases of new 
species. These changes are likely attributable to the over- 
all large increases in shrub and graminoid abundance 
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Table 4. The mean overlapping cover (%) in 2004 and 2011 of the 10 most typical species within each sampling unit and level of replication that 
contribute to the differences seen in the ordinations and ANOSIM analyses and the contribution (%) of each species to the total species pool in 
each group, using the SIMPER routine (Primer ver. 6). 



Upper SAS/summits 



Lower SAS/summits 



Species 


Family 


Life form 


2004 Mean 
abundance (%) 


2011 Mean 
abundance (%) 


Contril 


Kunzea muelleri 


Myrtaceae 


Shrub 


3.08 


3.44 


9.38 


Epacris microphylla 


Ericaceae 


Shrub 


3.76 


4.24 


9.03 


Poa sp. 


Poaceae 


Graminoid 


7.4 


9.76 


8.33 


Celmisia costiniana 


Asteraceae 


Forb 


3.72 


3.8 


5.69 


Empodisma minus 


Restionaceae 


Graminoid 


1.6 


2 


4.88 


Phebalium ovatifolium 


Rutaceae 


Shrub 


1.6 


1.52 


4.6 


Grevillea australis 


Proteaceae 


Shrub 


1.8 


1.48 


4.59 


Trisetum spicatum 


Poaceae 


Graminoid 


2.04 


1.88 


4.45 


Prostanthera cuneata 


Lamiaceae 


Shrub 


1.32 


1.32 


3.85 


Microseris lanceolata 


Asteraceae 


Forb 


1.44 


1.44 


3.75 


Kunzea muelleri 


Myrtaceae 


Shrub 


5.24 


5.72 


7.13 


Phebalium ovatifolium 


Rutaceae 


Shrub 


3.8 


3.72 


4.88 


Grevillea australis 


Proteaceae 


Shrub 


4.28 


3.48 


4.8 


Empodisma minus 


Restionaceae 


Graminoid 


4.6 


5.44 


4.4 


Prostanthera cuneata 


Lamiaceae 


Shrub 


3.2 


3.56 


4.22 


Epacris microphylla 


Ericaceae 


Shrub 


8.32 


9.68 


4.17 


Pentachondra pumila 


Ericaceae 


Shrub 


3.16 


2.92 


4.08 


Agrostis sp. 


Poaceae 


Graminoid 


0.04 


3.8 


4.05 


Carex breviculmis 


Cyperaceae 


Graminoid 


0 


3.84 


4.02 


Poa sp. 


Poaceae 


Graminoid 


12.64 


15.68 


3.31 


Kunzea muelleri 


Myrtaceae 


Shrub 


2.32 


2.72 


7.85 


Poa sp. 


Poaceae 


Graminoid 


7.6 


9.92 


7.6 


Celmisia costiniana 


Asteraceae 


Forb 


4.76 


4.6 


6.72 


Epacris microphylla 


Ericaceae 


Shrub 


3.52 


4.08 


6.24 


Phebalium ovatifolium 


Rutaceae 


Shrub 


2.04 


1.8 


5.57 


Empodisma minus 


Restionaceae 


Graminoid 


1.76 


1.8 


5.01 


Grevillea australis 


Proteaceae 


Shrub 


1.92 


1.6 


4.5 


Trisetum spicatum 


Poaceae 


Graminoid 


2.12 


2.2 


4.31 


Pentachondra pumila 


Ericaceae 


Shrub 


1.4 


1.6 


3.91 


Prostanthera cuneata 


Lamiaceae 


Shrub 


1.12 


0.96 


3.04 


Kunzea muelleri 


Myrtaceae 


Shrub 


4.72 


5.28 


5.82 


Phebalium ovatifolium 


Rutaceae 


Shrub 


4.52 


3.96 


5.25 


Empodisma minus 


Restionaceae 


Graminoid 


4.72 


4.32 


4.62 


Grevillea australis 


Proteaceae 


Shrub 


4.56 


3.72 


4.21 


Agrostis sp. 


Poaceae 


Graminoid 


0 


3.64 


3.76 


Carex breviculmis 


Cyperaceae 


Graminoid 


0.12 


3.72 


3.59 



Continued 
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Table 4. Continued 



Sampling unit/replication Species 


Family 


Life form 


2004 Mean 
abundance (%) 


2011 Mean 
abundance (%) 


Contribution (%) 


Pentachondra pumila 


Ericaceae 


Shrub 


3.92 


4.4 


3.51 


Epacris microphylla 


Ericaceae 


Shrub 


7.72 


9.04 


3.02 


Luzula alpestris 


Juncaceae 


Graminoid 


0.04 


2.92 


2.89 


Acetosella vulgaris 


Polygonaceae 


Forb 


1.16 


3.56 


2.87 



across the gradient, which in turn has consequences for 
driving overall community assembly processes given the 
trait-space and life-history strategies of these broad life- 
form groups (Diaz and Cabido 2001; Westoby et al. 2002). 
Although temperatures did not change significantly over 
the study period, precipitation did. The recent spikes in 
precipitation, however, are unlikely to have caused the 
patterns observed in the functional composition in this 
study; instead these may be realized in the subsequent 
1-2 growing seasons as the herbaceous vegetation re- 
sponds to the additional rainfall after many years of 
drought conditions. 

Community composition and ecology across 
the environmental gradient 

Our analysis of functional patterns from leaf traits may be 
used to help describe the mechanisms driving commu- 
nity composition across this elevation gradient of sum- 
mits (McGill et al. 2006). Overall, species have sorted 
themselves preferentially across the gradient, most likely 
as a result of abiotic filtering (Venn et al. 2011), the harsh 
conditions at the higher summits compared with the 
lower summits, and possibly through the additional impli- 
cations of biotic interactions. Both interaction types are 
ultimately predetermined by species' inherent functional 
traits and the functional trait space that each species oc- 
cupies (Diaz and Cabido 2001). Whereas a significant 
change in the CTWM across the gradient for a specific 
trait could indicate a threshold for an abiotic filter, it 
could also be due in part to biotic filtering, where positive 
or negative interactions between plants interact with en- 
vironmental factors preferentially across a gradient to ei- 
ther constrain or facilitate species with certain traits or 
trait combinations (Schob et al. 2012; Spasojevic and 
Suding 2012). 

Although we did not formally address the issue of intra- 
specific trait variation in this study, we noted that there is 
little variation in the traits measured within species 
across the gradient of sites; rather most of the variation 
in plant traits across the gradient comes from the vari- 
ation in species abundances, as demonstrated by the 



CTWM. Across gradients such as this one, there is an ex- 
pectation that the CTWM of each trait will be a reflection 
of a gradient of abiotic stress: higher at the higher- 
elevation summits and lower at the lower-elevation sum- 
mits. Conversely, the CTWM for certain traits might reflect 
some overriding biotic interactions which ameliorate 
stressful abiotic conditions as per the stress-gradient 
hypothesis (Bertness and Callaway 1994; Kikvidze et al. 
2006): in areas where abiotic stress is low and negative 
biotic interactions (competition) are expected (at the 
lower-elevation summits: higher temperatures, longer 
growing seasons) (Callaway et al. 2002; Cornwell and 
Ackerly 2009; Venn et al. 2009). This may explain the pat- 
terns in two of the four traits considered (the high CTWM 
values for SLA and leaf area), indicating that biotic filter- 
ing likely due to facilitation may be operating at the 
higher-elevation summits, allowing for relatively high 
productivity. Further experimental work in situ including 
measures of below-ground plant interactions will be ne- 
cessary to unravel the interacting effects of biotic and 
abiotic filtering (Wilson 1993) at these sites. At the 
lower summits, we observed relatively higher values in 
the CTWM for plant height and LDMC, strongly driven by 
the high proportion of shrub and graminoid species. 
A high proportion of these species can indicate fast and 
high biomass accumulation, lateral spread and competi- 
tive superiority (Westoby et al. 2002; Cornelissen et al. 
2003), as well as leading to increases in soil organic 
matter and available nitrogen (Myers-Smith et al. 2011; 
Vankoughnett and Grogan 2013). Here, plants may also 
be less constrained by (relatively) harsh conditions, but 
biotic filtering such as competition may override other 
assembly processes. 

Trade-offs between an energy-efficient strategy for 
nutrient conservation, with a high capacity for resource 
acquisition under harsh conditions, will likely drive future 
shifts in functional traits and species composition across 
the gradient (Choler 2005). However, the interaction be- 
tween abiotic filters (such as growing season length 
and variability) and biotic filters will ultimately determine 
whether short-term temporal shifts in species 
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Figure 3. The CTWM + 95 % confidence intervals for the four functional traits, plant height, leaf area, LDMC and SLA, based on species abun- 
dance in the upper and lower SASs. 2004 data = closed circles, 2011 data = open circles. Simple linear relationships (using all site aspect data) 
between site elevation and CTWM values are only displayed for R 2 > 0.5 and P < 0.05. 
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composition will eventually manifest as significant 
changes to observational trait patterns and ecosystem 
functioning in the long term. In the short term, the 



herbaceous species with high SLA and high productivity 
at high elevations must continue to 'hold their ground'. 
However, the relatively tall graminoids and shrubs, with 
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thick, tough leaves (low SLA) that dominate the lower 
summit areas and at lower elevations more generally, 
are likely to maintain their dominance, increase in 
abundance and out-compete other life forms in the 
local species pool. 
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